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The leading-order electron-atom bremsstrahlung is investigated within the rigorous relativistic 
approach based on the partial-wave representation of the Dirac wave functions in the external 
atomic field. Approximating the atomic target by an effective local potential, we calculate the 
double-differential cross section and the polarization correlations in a wide range of the impact 
energies. Connection between the bremsstrahlung at the hard-photon end point of the spectrum 
and the continuum-threshold limit of the radiative recombination is studied. A detailed analysis of 
the screening effect and the energy dependence of the polarization correlations is presented, with 
the main focus on the high impact energy region. 
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I. INTRODUCTION 



Bremsstrahlung, the emission of a photon by an elec- 
tron scattering from an atom, is one of the fundamen- 
tal processes that occur in the electron-atom collisions. 
For a large region of the impact energies of the incoming 
electron, the atom can be well represented by a static 
central (screening) potential, thus ignoring the exchange 
interaction between the incoming electron and the target 
and virtual excitations of the target. The corresponding 
approximation is sometimes referred as the "ordinary" 
bremsstrahlung mechanism, to distinguish from the so- 
called "polarization" bremsstrahlung, in which the ex- 
cess energy is transferred to the target and the pho- 
ton is emitted from the atomic core. The polarization 
bremsstrahlung is the dominant mechanism in, e.g., the 
proton-atom collisions (see, e.g., Ref. [l|), where the ordi- 
nary bremsstrahlung is suppressed by the small electron- 
to-proton mass ratio. In the electron-atom collisions, 
however, the polarization mechanism is usually not sig- 
nificant for the not-too-small energies of the incoming 
electron. 

Within the screening-potential approximation, the 
bremsstrahlung process can be described rigorously (to 
the leading order in the fine-structure constant a) within 
the exact relativistic approach based on the partial- wave 
representation of the Dirac continuum electron states 
with a fixed value of the asymptotic momentum. This 
approach is clearly preferable as compared to numerous 
approximate treatments reported in the literature but it 
is also more difficult for practical implementations. The 
partial-wave decomposition of the initial and final elec- 
tron states, together with the multipole expansion of the 
wave function of the emitted photon, lead to a large num- 
ber of expansion terms that have to be summed up until 
convergence is reached. Despite technical difficulties en- 
countered, the first partial-wave expansion calculations 
were reported already in 1960s [2, 3]. The first accurate 



numerical results were obtained by Tseng and Pratt in 
1970s. In a series of calculations 0t3|i they reported re- 
sults for the single- and double-differential cross sections 
for a wide interval of the impact energies. Polarization 
correlations in the double-differential cross section were 
studied in Ref. Q. More recently, several calculations of 
the triple-differential cross section and the corresponding 
polarization correlations were performed f9T-fllj|. 

Experimental investigations of the electron-atom 
bremsstrahlung, numerous in 1960s and 1970s [HJ], be- 
came rather sparse during the following decades. In the 
last years, however, advent of a new generation of the 
Compton polarimeters [l3[ revived the experimental in- 
terest to bremsstrahlung. In particular, it became possi- 
ble to perform experiments with the spin-polarized elec- 
trons, which opens up new vistas for experimental studies 
of various correlations between the polarizations of the 
incoming electron and the emitted photon. Such exper- 
iments are currently underway in Gesellschaft fur Schw- 
erionenforschung (GSI) [14j and Technical University of 
Darmstadt [l5j . 

Because of the experimental limitations in the past, 
the previous bremsstrahlung calculations were focused 
mainly on the cross sections. The only detailed partial- 
wave expansion study of the polarization correlations was 
accomplished by Tseng and Pratt [8|. Detailed as it is, 
this calculation is not sufficient to cover all present ex- 
perimental needs. In particular, the region of the impact 
energies E > 500 keV, presently available for the ex- 
perimental investigation in GSI, has not been carefully 
investigated in that work. Besides this, the results re- 
ported in Ref. |8| have never been carefully checked by 
an independent calculation. In the present investigation, 
we aim to cover this gap by checking and extending the 
previous bremsstrahlung results by Tseng and Pratt. 

The paper is organized as follows. In Sec. [Til wc 
present relativistic formulas for the double-differential 
bremsstrahlung cross section and polarization correla- 
tions, obtained within the density-matrix formalism. 
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Sec. IIIII describes the numerical approach used in this 
work. In Sec. IIVI we discuss the connection between the 
bremsstrahlung at the hard-photon end point of the spec- 
trum and the continuum-threshold limit of the radiative 
recombination. Numerical results are presented and dis- 
cussed in Sec. |V] 

The relativistic units (h = c = 1) are used throughout 
this paper. 



II. THEORY 

In the present investigation we consider the 
bremsstrahlung of an electron scattering from an 
atom, which is represented by a static central potential. 
The scattered electron in the final state is assumed to be 
not observed. The kinematics of the process is defined 
as follows. The reference frame is the rest frame of 
the atom. The z axis is directed along the asymptotic 
momentum of the incident electron pi. The xz plane 
(also referred to as the reaction plane) is defined by Pi 
and the momentum of emitted photon k. In this frame, 
the direction of the emitted photon is defined just by 
the polar angle 9k, cos 9 k — pi ■ k, where x = x/\x\. 

In order to study the polarization correlations between 
the incident electron and the emitted photon, we express 
the density matrix of the final state ( x k\\pf\kX'') 1 in terms 
of the initial-state density matrix as 

(fcA|p/|fcA') = ^2 J ' dn f (p l m l \a-uxe tk ' r \p f mf)* 

K) 
(i) 



x (pim^a ■ uy e^lpfuif) (pim^p^Pim'i) 



where Pi, rrii (rnQ and pf, rrif are the asymptotic mo- 
mentum and the spin projection of the incident and 
scattered electron states, respectively; k and A (A') are 
the momentum and the helicity of the emitted photon 
(A = ±1), u\ is the unit polarization vector of the pho- 
ton, and fi/ is the solid angle of the scattered electron. 
The energy of the emitted photon k = \k\ is fixed by 
k = Ei — ef , where e„ = ^Jp\ + m 2 . 

The advantage of the present formulation is that the 
density matrix (fTJ) contains the full information about 
the polarization properties of the emitted photon (see, 
e.g., Ref. 0), 
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(kX\p f \kX') = - Tr[p f ] 



1- 
Pi 



P 3 Pi - iP2 
-iP 2 1 - P 3 



(2) 



where P t are the Stokes parameters. We note that the 
sign of the Stokes parameters depends on the definition of 



1 More exactly, this is the reduced density matrix of the emitted 
photons, with the quantum numbers of the (unobserved) final- 
state electron traced out. In the context of the present paper, 
we refer to it just as the density matrix of the final state. 



the circular polarization unit vectors, see Eq. (fl4|) below 
and the text after it. The trace of the density matrix 
is, up to a prefactor, the double-differential cross section 
summed over all photon polarizations. Following Ref. [H, 
we introduce the normalized cross section a(k) as 



k da 

^ k ^ziTk= 27T 



d(cos ( 



da 



Z 2 dk dQk 



A-^£ d( cos, fc )Tr N , 



(3) 



where Z is the nuclear charge of the target atom and 
the continuum electron wave function is assumed to be 
normalized on the energy scale. 

The Stokes parameters Pi and P2 describe the linear 
polarization of the emitted photon. They can be deter- 
mined experimentally by measuring the intensities 1$ of 
the linearly polarized photon emission at different angles 
<f> with respect to the reaction plane, 



Pi 

Pi 



Jno — I f 



Ml* 



Iq° + 190 
/450 — /1350 



'45° 



:!5° 



(4) 
(5) 



Instead of the Pi and P2 parameters, it is often conve- 
nient to describe the linear polarization by the polariza- 
tion ellipse in the plane perpendicular to the photon mo- 
mentum k. The ellipse is defined by the relative length 
of the principal axis Pl which reflects the degree of the 
linear polarization, 



P L = X IP( + P$, 



and the tilt angle x. 



X^acrtan^ 



(6) 



(7) 



The Stokes parameter P3 defines the degree of the circu- 
lar polarization of the emitted photon, 



Pi 



W(+l) - W(-l) 
W(+l) + W(-l) ' 



(8) 



where W(+l) and W(—l) are weights of the right- and 
left-polarized photons, respectively. 

We now turn to the evaluation of Eq. (fTJ). The initial- 
state density matrix is represented in terms of the spher- 
ical tensors p£ of rank k = and 1 as [l6| 

<P i m i |p i |pX> = E(- 1 ) V2 " m ^vL i ,V 2 -^^- ( 9 ) 



The components of p$ are expressed in terms of 
the polarization vector of the incoming electron Pi — 
(Pr, Py, Pz) as 

(10) 
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The spherical- wave decomposition of the ingoing (+) and 
outgoing (— ) continuum electron wave function with a 
definite asymptotic momentum is 



\pm) = 4ir 



/ 1 € 



where j = \k\ - 1/2, I = \k + 1/2| - 1/2, A K = 
&k + (I + 1) t/2, and tr K is the asymptotic phase of the 
wave function, see Eq. (fTBl below. The function \eKfi) is 

the Dirac eigenstate with the energy e = y'p 2 + m 2 , the 
rclativistic angular quantum number re, and the angular 
momentum projection /i, represented by 



\fku\-1 9eA r )x^(r) 
l£K ^-\ if e , K (r) X -^(r) 



(12) 



where g and / are the upper and lower radial compo- 
nents, respectively, and x K( i are the spherical spinors [18| . 
The wave function is normalized on the energy scale and 
its asymptotic (as r — > oo) behavior is 



, s 1 ( e + m 

geA r ) - - 

r \ irp 



1/2 



cos \pr + a K + rj ln(2pr)l , 

(13) 



where r) = Z^ae/p, and Z^ is the effective nu- 
clear charge of the atom at large distances, Zoo = 
lmv-^oo rV SCI (r), with V scr being the atomic potential. 



The spherical-wave decomposition of the photon held 
with the helicity (circular polarization) A = ±1 is 

u x e ikr = V2^J2 * L V2LTl(i\r a ( [ ] M (r) D L MX (k) , 

LMp 

(14) 

where the components of the circular-polarization vector 
U\ are defined as [l8| u\ = (u x + iu y )/y/2 and u_i = 

(u x — iUy)/V%, -Dma i s Wigner's D function [ijj], a^ M 
are the magnetic (p = 0) and electric (p = 1) vectors 
defined by 



a( LM{r) = jL(kr)Y LLM (r) , 



(15) 



"ImW =3L-i(kr) 



j L+1 {kr) 



L + l 
2L + 1 



- L L — l M 



(f) 



L 



■LL+1M 



(r) , (16) 



2L + 1 

Yjlm 1S the vector spherical harmonics [l9j . 

Y JLM (r) = C J Z,Xa y Lm (r) e a , (17) 

ma 

and e a is the spherical component of the unity vector. 

Inserting expansions ©, CU} and flU into Eq. ((TJ and 
employing the standard angular-momentum technique, 
we arrive at 



(k\\ Pf \k\') = 8(2^ y, E E^(fc)pSU i,1_,i " i+L ' eiA ^"^ [^^j*.^'*. 1 *'^": 



1/2 



KiK f Kf LL' Kgt 7172 



^ ^ U L'A',L-A u iiO,ZJO Ly s-7i,«7i 1 j' g Jj 



i/a J/a re 

ii ii 5 



pp' 



a • a 



(p) 



c", ti 



(p') 



£ f K f 



(18) 



where [xi, x 2 , ■ ■ ■} = (2a;i + l)(2x 2 + 1) . . ., l a = |« a + 

V 2 1 - V 2 ' C = l K i - V 2 I - y 2 ^ ia = |«a| - x /2, and the 
reduced matrix elements are evaluated in Appendix [X] 

In this work, we will present our results in terms 
of the differential cross section and the Stokes param- 
eters of the emitted photon as functions of the polar- 
ization vector of the incident electron: <r(P x , P y , P z ), 
P!(P X ,P V ,P Z ), P 2 {P x ,P y ,P z ), and P 3 {P x ,P y ,P z ). We 
consider four choices of the polarization of the inci- 
dent electron: (i) unpolarized, (P x , P y , P z ) = (0,0,0); 
(ii) polarized transversely within the reaction plane, 
(P x ,P y ,P z ) = (1,0,0); (iii) polarized perpendicularly to 
the reaction plane (P x , P y , P z ) — (0,1,0); and (iv) lon- 
gitudinally polarized, (P x ,P y ,P z ) = (0,0,1). Not all of 



the polarization correlations for these 4 choices are in- 
dependent and nonvanishing. Properties of the angular- 
momentum coefficients in Eq. (TT8"]) lead to the following 
identities, 

Pi (0,0,0) = Pi (1,0,0) = Pi (0,0,1), (19) 
P 2 (0,0,0) = P 2 (0,1,0) = 0, (20) 
P 3 (0 ) 0,0) = P 3 (0 ) 1,0) = 0, (21) 
da(0, 0, 0) = da(l, 0, 0) = da(0, 0, 1) , (22) 

where da(P x , P y , P z ) = (k/Z 2 ) d<r(P x , P y , P z )/ {dk rffi fc ). 
Taking into account the above identities leaves us with 
8 independent quantities to be calculated: da(0, 0,0), 
(2(7(0,1,0), Pi (0,0,0), Pi (0,1,0), P 2 (1,0,0), P 2 (0,0,1), 
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P 3 (1,0,0), andP 3 (0,0,l). 

In previous calculations reported in the literature, 
the polarization correlations were often parameterized in 
terms of the coefficients Cy introduced by Tseng and 
Pratt [8] . In order to simplify comparison with previous 
studies, we give the list of correspondence between the 
coefficients Cy and the present notations, 



C 03 = Pi (0,0,0), 

Cxi = -P 2 (1,0,0), 

C 12 = -P 3 (1,0,0), 

C 23 = Pi (0,0,0) -Pi (0,1,0). 

C 3 i = P 2 (0,0,1), 

c 32 = P 3 (0,0,1), 

d<x(0,l,0) 
20 da(0,0,0)' 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 

(29) 



III. NUMERICAL CALCULATION 

The problem of calculating the final-state density ma- 
trix (|18p is now reduced to the evaluation of the radial in- 
tegrals and the summation over the angular-momentum 
and multipole quantum numbers. The radial integrals to 
be evaluated are 



J," (a, 6) 



drr , g ea ^.{r)f Sb ^{r)ji{kr), (30) 



Jf L (a,b) = / drr'f eailta (r)g eb , Ki (r)ji(kr), (31) 
Jo 

where g and / are the radial components of the contin- 
uum Dirac wave function and ji is the spherical Bessel 
function. A straightforward numerical evaluation of these 
radial integrals is problematic since all three functions in 
the integrand are highly oscillating and slowly decreasing 
for large values of r. 

Several methods were reported in the literature for 
evaluation of such highly oscillatory integrals. For the 
case of the point-Coulomb potential, the integrals can 
be evaluated analytically in terms of the Appel functions 
of three complex variables @, [2(| ■ Weak points of this 
method are, first, the restriction to the particular choice 
of the atomic potential and, second, the absence of reli- 
able numerical methods for evaluation of the Appel func- 
tions for a wide range of parameters and arguments. 

The method used in calculations by Tseng and Pratt 
[3, H[ is based on dividing the integration region (0, oo) 
into two parts, the inner and the outer ones, in such a 
way that the wave function in the outer region can be 
approximated by the (phase-shifted) free-field solutions. 
In the inner part, the Dirac equation is solved and the ra- 
dial integrations are performed numerically. In the outer 
region, the analytical form of the free-field solution is 
exploited and the radial integrals are performed by the 
so-called integration by parts method, which was never 
described in the literature [2l|. 



In the present work, the radial integrals are evaluated 
numerically after rotating the integration contour in the 
complex r plane. For the first time this elegant method 
was used probably in Ref. [22j . in connection with nu- 
clear collision problems. For the case of the integrals 
involving three spherical Bessel functions, this method 
was studied in detail in Ref. [23]. Integrals involving the 
nonrelativistic Coulomb functions and negative powers 
of the radial coordinate were investigated in Ref. [24] . 
More recenly, the complex-plane rotation method was 
used in calculations of the (e, 2e) process [25| and the 
Compton scattering (26|. In the present work, we ex- 
tend the complex-plane integration method to the case 
of the Dirac Coulomb wave functions. The calculational 
scheme is described in detail in Appendices [B] and [C] 
This method allowed us to evaluate the radial integrals 
required up to typically 9-digit precision. 

The continuum-state wave functions for the general 
case of the screening potential were obtained by nu- 
merical solution of the Dirac equation with help of the 
RADIAL package by Salvat et al. [27| . For the point- 
Coulomb potential, we used the analytical representation 
for the Dirac eigenstates in terms of the Whittaker func- 
tions, see Eqs. (|B2j and (|B3]) . The Whittaker M and W 
functions were evaluated by a generalization of our codes 
developed in Ref. [Hj|. 

In the present work, we perform calculations for two 
choices of the atomic potential: (i) the point-nucleus 
Coulomb potential and (ii) the screening potential of a 
neutral atom. The screening potential was constructed 
as 



V SCT {r) = V nuc (r) + a 



dr' 



1 



max(r, r 



TTP(^), (32) 



where V nuc is the potential of the (extended-size) nucleus, 
Rq is the radius of the atom, and p(r) is the radial elec- 
tron density of atomic orbitals normalized by 



Rd 



dr p(r) = Z . 



(33) 



The radial atomic electron density was calculated by 
the multiconfigurational Dirac- Fock method by using the 
GRASP package [29[ . We note that the previous calcula- 
tions by Tseng and Pratt d, [|| used a somewhat simpler 
model of the screening potential, which included only a 
local part of the exchange interaction between the atomic 
electrons. 

Once the radial integrals are successfully evaluated, the 
main problem in numerical calculation is to ensure the 
convergence of the multiple sum in Eq. (|18[) . When all tri- 
angular selection rules are taken into account, two sum- 
mations remain unbound, which can be chosen to be |/tj| 
and | Kf | . The convergence properties of the sum depend 
strongly both on the kinetic energy of the incoming elec- 
tron E and on the fractional part of it carried away by the 
photon, k/E. The convergence is fast for E < 1 MeV and 
close to the hard photon end of the spectrum (the tip re- 
gion, k/E ~ 1) and becomes rather slow for larger impact 
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energies and(or) near the soft photon end. In the tip re- 
gion, the expansion over the final-state quantum number 
| At/ 1 converges rapidly, so that a non-symmetric partial- 
wave cutoff is appropriate. Our calculations in this region 
included up to (nj,n/) = (80,10) partial waves. Further 
away from the tip region, we used the symmetric con- 
figurations up to (rij,n/) = (55,55) partial waves, with 
further extension being complicated by numerical insta- 
bilities. 



IV. CONNECTION BETWEEN 
BREMSSTRAHLUNG AND RADIATIVE 
RECOMBINATION 

From the general physical point of view, it is natural to 
expect a relation between the bremsstrahlung (BS) at the 
hard photon end of the spectrum (the tip region) , where 
the electron transfers all its kinetic energy to the photon, 
and the radiative recombination (RR) into highly excited 
atomic states. This connection was studied previously by 
a number of authors (30l - [33l [ . Approaching the tip region 
from the continuum side was invesigated in Ref. [3~j| . In 
this work, we would like to demonstrate the connection 
between these two processes by explicit numerical calcu- 
lations. 

For the Stokes parameters, the connection between RR 
and BS is most transparent. It can be observed that, for 
RR into atomic states with a given angular momentum 
quantum number At, the Stokes parameters have a well- 
defined limit as the principal quantum number n — > oo. 
In order to make a connection to BS, one has to account 
for the recombination into subshells with different values 
of At (for a fixed value of n) and then evaluate the limit 
of n — > oo. More specifically, the RR Stokes parameter 
Pi is extrapolated to the continuum threshold as 



ptip 



lim 



J2 K [h°{nn) - J 90 o (nn)} 
Y, K [ J o»H +^90° (n«)] 



(34) 



where I x (n/c) denotes the intensity of the photons emit- 
ted under angle x m RR °f an electron into a bound 
atomic substate |nAt). For explicit formulas and details 
on I x {nn), we refer the reader to the recent review by 
Eichler and Stohlker [l7|. The other Stokes parameters 
Pi and P3 are obtained analogously to Eq. (jM)) , by using 
Eqs. © and ©. 

In actual calculations, the limit n — > 00 was evaluated 
by going to sufficiently high values of n (typically, n ss 
15 — 20), whereas the summation over k was restricted 
to a few terms. (The s, p, d, /, and g waves were taken 
into account.) 

A somewhat different procedure is required for ap- 
proaching the continuum threshold limit in the case of 
the cross section. It is well-known that the cross sec- 
tion of RR into the individual Rydberg states falls off as 
r~> n~ 3 . Making a connection to the BS 33], one should 
take into account that the individual levels cannot be re- 
solved experimentally as n — > 00. The measured quantity 



depends on the resolution of the photon detector Afc and 
is given by 



\ dn k / Ak Ak ^ dSlk 



(35) 



n>no 



where cr^? is the cross section of RR into an individual 
level with given n and k and uq is defined by the condition 
that the binding energy E no = Ak. In order to approach 
the continuum threshold, we sum over At and take the 
limit uq — > 00, 



dk dQi. 



lim 



n>no K J 



Ak(n ) 



(36) 



In the limit n — > 00, one can use the nonrelativistic ex- 
pression for the binding energy, Ak(n) — m(Za) 2 /(2n 2 ), 
and the well-known asymptotic behavior of the RR cross 
section, da^/dilk = A K (6)/n 3 , with the result 



dk dQi- 



y 

^-^ m(Za) 2 



(37) 



In our calculations, we determine the parameter A K (6) 
by fitting the RR cross sections calculated for a series of 
n, to the 1/n 3 scaling law. 

Fig. Q] illustrates the connection between BS and RR, 
as obtained in our numerical calculations for the dou- 
ble differential cross section and the Stokes parameters 
in the case of the point Coulomb potential. We observe 
good agreement in all cases; the small remaining devia- 
tion is attributed to extrapolation errors. It is important 
that our RR and BS calculations are completely indepen- 
dent from each other. The BS calculation is performed 
as described in this work, whereas the RR calculation 
is carried out as reported in Ref. [35j ]. with the radial 
integrations performed analytically. So, the agreement 
observed is also as an important cross-check between two 
different calculational approaches. 



V. RESULTS AND DISCUSSION 

We begin this section by comparing our numerical re- 
sults with those obtained in the best previous partial- 
wave expansion calculations by Tseng and Pratt [U, H[ . 
TableUpresents comparison for the bremsstrahlung cross 
section a(k) given by Eq. §5§ for the gold target, both for 
the Coulomb and screening potentials. We observe ex- 
cellent agreement between the calculations performed for 
the Coulomb potential. For the screening potential, all 
our results are by about 1% smaller than the ones by 
Tseng and Pratt. This difference is apparently due to a 
more realistic screening potential used in this work [see 
the discussion after Eq. ([33])]. Results of the two calcula- 
tions for the double differential cross section and for all 
independent polarization corrections Cy [see Eqs. (|23[1 - 
|] are compared in Fig. [5] for the case of the Coulomb 
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potential. We find nearly perfect agreement in all cases. 
As a matter of fact, the agreement observed is quite re- 
markable, taking into account that the calculation by 
Tseng and Pratt was accomplished four decades ago. 

In Fig. [3J we study the effect of the screening on the 
cross section and the Stokes parameter Pi for the ini- 
tially unpolarized electrons. First of all, we observe that 
the sign of the effect varies: the screening reduces the 
cross section but increases the Stokes parameters. For 
the cross section, the screening effect depends strongly 
on the nuclear charge Z . The effect is barely recogniz- 
able (in the region of energies studied) for carbon, while 
for gold, it is significant for all energies. 

It is interesting to observe that the bremsstrahlung on 
the bare nucleus and the neutral atom are very much 
alike for not-too-small impact energies. This is explained 
by the fact that the dominant contribution to the ma- 
trix element of the amplitude of the process comes from 
the small distances of configuration space. At these dis- 
tances, the continuum-state electron "feels" mainly the 
bare nuclear charge and the screening effect is relatively 
small. At larger distances, the continuum-state electron 
wave function oscillates rapidly, so that the contribution 
to the radial integal is small. When the impact energy 
increases, the region of configuration space responsible 
for the dominant part of the integral shifts towards the 
nucleus and the screening effect decreases. 

For E = 20 kcV and gold target, the screening reduces 
the forward-scattering cross section by more than a fac- 
tor of 2, whereas at energies of a few MeV, it is still a 
10% effect. For the Stokes parameters, the screening ef- 
fect is smaller than for the cross section and stays within 
10% even at the lowest energy studied, E = 20 keV. 
We observe that for energies below 20 keV and gold tar- 
get, calculations of the cross section within the screening- 
potential approximation rapidly become meaningless as 
the screening effect tends to grow very fast as the energy 
decreases. On the contrary, the Stokes parameters turn 
out to be much less sensitive to the electronic structure of 
the target. We also find that for the high energies above 
1 MeV, the screening effect is significant mainly for the 
cross section at close-to-zero angles, where it amounts to 
about 10%. 

It is interesting to observe that the dependence of the 
forward-scattering cross section on the impact energy in 
the tip region is drastically different for the \ow-Z (car- 
bon) and high-Z (gold) targets. This difference is studied 
in more detail in Fig. 2] We observe that for carbon, the 
forward-scattering cross section at the tip nearly vanishes 
for large range of impact energies, as could be expected 
from the nonrelativistic theory. On the contrary, the gold 
target corresponds to the highly relativistic regime and 
its forward-scattering cross section is strongly non-zero 
for all energies. At high impact energies E > 0.5 MeV, 
the forward-scattering cross section grows rapidly with 
the increase of E, whereas the angular distrubution of 
the cross section becomes localized at the increasingly 
smaller regions near the forward direction. 



In Fig. [5l we study the energy dependence of the 
Stokes parameter Pi, the only nonvanishing Stokes pa- 
rameter for the initially unpolarized electrons. We ob- 
serve that the shape of the angular distribution of Pi de- 
pends strongly on the energy of the initial electron E. For 
small E, Pi tends to its nonrelativistic value of Pi = 1, 
whereas for the increasingly larger energies, maximum of 
the angular distribution decreases and shifts towards the 
forward direction and its tail crosses the zero axis (the 
"crossover" feature @ ) . On the contrary, when the initial 
energy E is fixed and the energy of the final electron is 
varied, the shape of the angular distribution stays much 
the same but its amplitude decreases. 

In Fig. |6j we plot the Stokes parameter Pi as a func- 
tion of the fractional energy carried away by the photon, 
k/E, for two targets, carbon and gold, and for two values 
of the angle of the emitted photon. These plots corre- 
spond to the experiment that is currently underway in 
TU Darmstadt [T^. 

We now turn to the case of the initially polarized elec- 
trons. In the present investigation, we restrict ourself to 
studying the longitudinal polarization. Fig. [7] presents 
our results for the Stokes parameters P2 and P3 as func- 
tions of the initial energy E and of the energy of the final- 
state electron. (We recall that Pi for the longitudinally 
polarized electrons is the same as for the unpolarized 
ones and shown in Fig. [SJ) The Stokes parameter P2 is 
of pure relativistic origin and thus is very small for small 
initial energies. With increase of E, P2 becomes com- 
parable to Pi, has its maximum at about E = 1 MeV, 
and then gradually descreases. The third Stokes param- 
eter P3 is also a relativistic effect and thus vanishes for 
small initial energies. However, in the high-energy region 
and k/E ~ 1, P3 changes its behaviour drastically and 
approaches unity everywhere except for the backward di- 
rection. 



VI. SUMMARY 

In this paper, we report a detailed study of the 
electron-atom bremsstrahlung process within the rigor- 
ous relativistic approach based on the partial-wave ex- 
pansion of the Dirac wave functions in the external 
atomic field. Assuming that the final-state electron is 
not observed, we evaluate the double-differential cross 
section and all polarization correlations. Unlike in the 
previous studies, our description of the polarization cor- 
relations is formulated entirely in terms of the Stokes 
parameters, which are directly related to quantities ob- 
served in modern experiments. For our calculations, we 
developed an efficient and reliable scheme of evaluation 
of the radial integrals for the free-free transitions, based 
on the complex-plane rotation of the integration contour. 
The method is applied for the Dirac solutions in both the 
point-Coulomb potential and the finite-range screening 
potential. 

The numerical procedure was carefully checked by 
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comparing our results against those reported in the lit- 
erature. Comparison of our results obtained for the 
hard-photon end point of the bremsstrahlung spectrum 
with the extrapolation of the radiative recombination re- 
sults yielded a numerical proof of the connection between 
bremsstrahlung and radiative recombination and served 
as an additional check of our computational scheme. 

The numerical results reported present a detailed anal- 
ysis of (i) the screening effect induced by the electrons 
of the target on the cross section and polarization cor- 
relations and (ii) the energy dependence of the polar- 
ization correlations, with the main focus on the high- 
energy region, which is of primary interest in the future 
experiments at GSI. We conclude that the tip region of 



the bremsstrahlung spectrum is the most appropriate for 
studying the polarization correlations, as all polarization 
correlations have their maximum values there. 
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Appendix A: Reduced matrix elements 

Results for the reduced matrix elements in Eq. (p~8|) are 



(-i)(e a K a \ 



a ■ a 



(0) 



e b K b ) = J ; 12 (a, b) s u (K a ,-Ki) - J; 21 (a, 6) s u (-K a ,K b ) . 



(Al) 



{-i)(e a K a \\a ■ a^H^Kb) = J ^j^j 



21 + 1 



Jl^a.b) su-x(K a ,-Kb) - Ji\(a,b) su-x(-K a , K b ) 
Ji+i(a, b) sii + i(n a , -Kb) - Jfli{a, b) su + i(-K a , Kb) 



where the angular coefficients are given by 



(A2) 



r~£~ I ji h V 2 I 

fl£j(«H,Ka) = <Ki||o--ybj||/sa) = (ki\\[Yj ® <t}l\\k 2 ) = \j — (-1)' 2 [ji,j2,k,k,L] 1/2 Cf° 0Aa0 < h h V 2 > , 



(A3) 



Yjlm are the vector spherical harmonics defined by 
Eq. (IT7|) and the radial integrals J ; 12 (a, b) and J 21 (a, &) 



are defined by Eqs. J30J and (l31j) . 



Appendix B: Free-free integrals: Coulomb case 



We need the radial integrals of the form 



J?(ei,Kl,£2,K2) 



drr 2 f euKll i(r) / E2 ,« 2 ,j(r) ji(kr) , 

(Bl) 



where ji is the spherical Bessel function, k = \si — £2\, 
and / EiKjl ee g e K and / e:Ki2 = / e , K are the upper and the 
lower radial components of the Dirac wave function. In 
the following, we will assume that e\ > £2 > m, w hich 
entails that pi — p 2 — k > 0, where pi — \J £ 2 — rn 2 . By 
this restriction, we exclude the possibility £1 — £2, in 
which case the method described in this section is not 
applicable. 

For £ > m, the radial components of the Dirac- 
Coulomb wave functions normalized on the energy scale 
are given by [36| : 

g e , K (r) = AT K% /Fn(2pr)- 3 / 2 Re( e 4 ^-f(^)] 



feAr) 



x (7 + i v ) M_i / 2_ <I ,, 7 (2ipr)| , (B2) 
- N K V£~T (2pr ) - 3 / 2 Im j e l [ s « " * ( 7+ * )] 
x (7 + w?) Af_i /2 _ 4 , j7 (2ipr) j , (B3) 



where p = \J £ 2 — in 2 , 



Nk 



t] = Zct£/p, 7 = -\/k 2 — (Za) 2 , 



„2i<5„ 



P p ^/2 l r (7 

r(2 7 + i) 



iri/e 



7 + w; 



(B4) 



(B5) 



and M Qi/ 3 is the Whittaker function of the first kind [371 ]. 

It is clear that the integrals J\ 3 can be expressed in 
terms of the integrals involving two Whittaker functions 
and a Bessel function, 



l 0:1,71,02,72 



,i(Pi,P2,A;) 



dr r 1 M ai , 7l (2ipir) 



x M Q2:72 (2ip 2 r) j; [kr) , (B6) 



with the momenta satisfying the condition p\—p2~k > 0. 
The integrand is highly oscillatory and slowly decreasing 
function for large values of r, so a straightforward numer- 
ical evaluation of this integral up to a high accuracy is 
practically impossible. The method used in the present 
work [23j is based on the analytical continuation of the 
integrand into the complex r plane. More specifically, 
the integrand is separated into two parts which decrease 
exponentially in the upper or lower half of the complex 
plane. By appropriate rotations of the integration con- 
tour, the original oscillating integrand can be converted 
into two smoothly decreasing ones. The resulting inte- 
grals can be easily calculated by Gauss-Legendre quadra- 
tures up to typically 9 digit accuracy. 

To realize this algorithm, we represent the first Whit- 
taker function in the right-hand-side of Eq. (IB6|) in terms 
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e(deg) e(deg) e(deg) e(degJ 



FIG. 1: (Color online) Comparison of the bremsstrahlung at the hard photon end of the spectrum (solid line, black) with the 
continuum-threshold extrapolation of the radiative recombination (dashed line, red), for the double differential cross section 
(k/Z 2 ) da/(dkdQk) and the Stokes parameters. Calculations are performed for Z = 79, E — 100 keV, and the point-nucleus 
Coulomb potential. 




FIG. 2: Comparison of the present results (solid line) with those by Tseng and Pratt 4, 8] (open circles) for the double differential 
cross section (k/Z 2 ) da / (dk dQk) for the initially unpolarized electrons, in mbarn/sr, and for the polarization correlations CV, . 
The calculations are performed for Z = 79, E = 500 keV, k = 250 keV, and the point-nucleus Coulomb potential. 



of the Whittakcr functions of the second kind 1 37 



Bessel functions (for p and k > and \z\ — > oo) is 



r(2/3 + i) 



r(/3 + a + l/2) 



(B7) 



where s — 1 if Im(z) < and —1 otherwise. For r = 
R + iz, the asymptotic behavior of the Whittakcr and 



M ai/3 (2ipr) 



W a ,f}(2ipr) 



ji{kr) 



'(B8) 

where only the leading exponential behavior is kept. 
These results, together with the condition on the mo- 
menta, pi > P2 + k, show that the representation (|B7[) 
applied to the Whittakcr function with the largest mo- 
menta M aira (2ipir), splits the integrand of Eq. (|B6[) 
into two parts, one of which [with W- ai )71 (— 2ipir)] is 
regular in the upper half of the complex r plane and the 



10 



1 

b 



1 

C 





Au 


e = o u '/ 


. e = 6o 

\ 




/ 






// 
// 
// 

=sS$T 8 = 30" 




100 


1000 


\ ll 6() 




C 


\ \ 

e = 30" \\ 






^ e = 0" 








100 


1000 




E(keV) 







1 .0 




0.8- 




0.6- 











0" 


0.4- 
















0.2- 




0.0- 




-0.2- 




1.0 




0.8 - 




0.6 - 












0.4 












0.2- 




0.0- 




-0.2 








c 


//' e = 10" 


8 = 60" \ 




\ 8 = 30° 



E(keV) 



FIG. 3: (Color online) Double differential cross section (k/Z 2 ) da / (dk dtl^) and the Stokes parameter Pi for several angles of 
the emitted photon 9, as functions of the kinetic energy of the incoming electron E, for gold (upper row) and carbon (lower 
row) targets, for the Coulomb (dashed line) and screened (solid line) potentials. For all graphs, the kinetic energy of the final 
electron is Ef — 1 keV and initially unpolarized electrons are taken. 




6 (deg) e (deg) 



FIG. 4: (Color online) Double differential cross section (k/Z 2 ) da/(dk dO.^) for different values of the kinetic energy of the 
incoming electron E, for neutral carbon (left) and gold (right) targets. Kinetic energy of the final electron is Ef = 1 keV and 
initially unpolarized electrons are taken. 
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FIG. 5: (Color online) Stokes parameter Pi for the neutral gold target and initially unpolarized electrons. Left panel contains 
plots for different values of the kinetic energy of the incoming electron E and the kinetic energy of the final-state electron fixed 
by Ef = 1 keV. In right panel, the initial energy is fixed by E — 1 MeV and the final-state energy is varied. 




FIG. 6: Stokes parameter Pi for the initially unpolarized electrons with E — 100 keV, for different values of the energy k and 
the angle 9 of the emitted photon, for gold and carbon targets and for the Coulomb and screening potentials. 
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FIG. 7: (Color online) The same as Fig.[5]but for the Stokes parameters P2 (upper row) and P3 (lower row) and for the initially 
longitudinally polarized electrons. 



other [with W aini (2ipir)] is regular into the lower half- 
plane. In both cases, the resulting integrand falls off as 
~ e~( pi ~ P2 ~ fe 'l z l , which makes possible accurate numeri- 
cal evaluation of the corresponding integrals. 

In actual calculations, one should keep in mind that 
Eq. (IB7|) represents the regular (for small z) function 
M in terms of the irregular functions W. Because of 
this, it is advantageous to integrate Eq. (|B6p along the 
real axis up to a certain value of r = R and to perform 
the rotation of the integration contour for r > R. So, 
Eq. (|B6[) is represented as 



I = 



dr I(r) 



/>oo 

/ dz i [l + (R + iz) - I- (R - iz)] , 
Jo 

(B9) 



where J(r) stands for the integrand of the right-hand-side 
of Eq. (|B6I) . and I+(r) and I-(r) are the parts of I(r) 
regular in the upper and lower half-plane, respectively, 
/(r) = J + (r)+J_(r). 



Appendix C: Free-free integrals: neutral atom 

We now consider the potential in the Dirac equation 
to be a finite-range screening potential (having in mind 



a neutral atomic system), 

V(r) = V scr (r) , 



(CI) 



where V scr (r) — for r > Rq. In this case, in the outer 
region r > Rq, the solutions of the Dirac equation can be 
expressed in terms of the free Dirac cigenf unctions, given 
by 



ji(pr) cosS - yiipr) sin 5 



(C2) 



/e (0) (r) = \K\Jp{e-l) 



K V 7T 



jjipr) cos S — yj(pr) sin 5 



(C3) 



where ji and yi are the spherical Bessel and Neumann 
functions, respectively; I = \k + 1/2| — 1/2, I = \k — 
1/2 1 — 1/2, and 6 is the scattering phase shift induced 
by the finite-range potential V scl - and determined by the 
matching procedure at r = Rq. (In the absence of any 
field, 5 = 0.) Note that the definition of the overall sign 
in the wave functions agrees with that of Ref . [36| . The 
definition used in Ref. [27J differs from the present one 
by a factor of —k/\k\. 
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For the evaluation of radial integrals, it is convenient 
to express the free Dirac eigenfunctions in terms of the 
spherical Hankel functions of the first and second kind 
(/if and h\ f , respectively), 



fl$(r) 



p(e + l) 1 

7T 2 



h[ 1] {pr) e lS + hf ] [pr) e~ lS 



,(2) 



(CM) 
(C5) 



Taking into account that h\ (pr) 



if~\pr) ~ e lpr and hf'\pr) ~ 
e~ lpr as r — > oo, we conclude that the first terms in the 
brackets of Eqs. (IC4[) and (|C5[) are regular in the upper 
half of the complex r plane, whereas the second terms 
are regular in the lower half-plane. 



The integral involving the free Dirac eigenfunctions can 
be evaluated as 



J l l 1 (£ 1 ,Kl,S2,K2,R) 



dz 'i 



f(0) 



*2 AO), 



f(0) 



(C6) 



where x = R+iz, and superscripts "+" and "— " label the 
parts of the eigenfunction that are regular in the upper 
and lower half-plane, respectively. Using the property of 
the spherical Hankel functions 



hf\z) =hf\z*) 



(C7) 



the expression for J]j is reduced to 



dz x 2 



f 

JI j (ei,Ki,£ 2 ,k 2 ,R) =-2Im / 

Jo 

(C8) 

It can be easily seen that the integrand in the above 



expression falls off as e~( pi_P2 ~ fc " )z for large z. 

For r smaller than the radius of the atom Rq, the 
wave functions are obtained numerically by solving the 
Dirac equation, so this part of the radial integration has 
to be performed along the real axis. When the energy 
of the incoming electron becomes large, the integration 
within the radius of the atom (which is of order of sev- 
eral atomic units) might become troublesome as the in- 
tegrand is rapidly oscillating. To simplify the numerical 
evaluation, we exploit the fact that, for heavy atoms, 
the maximum of the screening potential V SCI is localized 
close to the nucleus. Further away from the maximum, 
the numerical solutions of the Dirac equation resemble 
the free asymptotic solutions, so that the difference be- 
tween them is a smooth, rapidly decreasing function. We 
thus represent the radial integrals JV [Eq. (|B1[) ] as 



J 



J/ J (ei,Ki,£ 2 ,K2) = / drr 2 f euKul (r) f e2 ^ 2 j(r)ji(kr) 



Ro 
IR 

- 21m 



dr r 



/ei, «!,<(?") .fs 2 -K,2,j{ r ) /ei,Ki,i( r ) fs2,K 2 ,j( T ) 



p(0) 



ji{kr) 



(C9) 



where x = R + iz. The free parameter R < Rq is chosen as the smallest distance for which the exact Dirac solution 
fe,K,,i( r ) resembles the free asymptotic solution j(r). 



